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Abstract 

Direct photon production in association with a heavy quark can serve as an excellent tool for the study of the heavy 
quark distributions. Currently it is believed that heavy quarks are produced radiatively inside the nucleus, and so there 
is no need to take into account heavy quark parameters inside global PDF analysis. Certain models taking into account 
the possibility of an intrinsic charm component exist. Here we present how these affect the y+c cross section. While at 
pA collisions the potential of this process to constrain the gluon nuclear PDF which is currently largely unconstrained 
is presented. 
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1. Introduction 

7.7. Direct y+ Q 

In order to gather more information on the heavy 
quark PDFs it is important to investigate in detail pro- 
cesses sensitive to them. Direct photon production in as- 
sociation with a heavy quark (c or b) is exactly one such 
process due to its dependence on subprocesses sensitive 
to the initial state heavy quark (HQ) distributions ||T]. 

At LO it has a simple form - it consists of one hard- 
scattering subprocess, g + Q ^ y + Q (Compton sub- 
process), and also fragmentation contribution^ these 
being greatly suppressed due to experimental isolation 
requirements. At NLO the number of contributing hard- 
scattering subprocesses increases to seven 

g + g^ Q + Q + y Q + Q^ Q + Q + y 

g + Q^g + Q + y Q + Q^ Q + Q + y 

Q + q^q+Q + y q + q^ Q + Q + y 

Q + q^ Q + q + y 

also including NLO fragmentation effects, which are 
small due to the isolation cuts as in the LO case. Given 
that most of the subprocesses are initiated by either 
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' The LO fragmentation contributions consist of all 2 — » 2 subpro- 
cesses of order <5(aj) containing at least one heavy quark in the final 
state convoluted with the photon fragmentation function Dyi^ g(z, Q^). 

^This calculation is performed in the variable flavor number 
scheme, where the heavy quarks are treated as massless. 



heavy quarks or gluons, and in the proton the gluon dis- 
tribution is well known, it is possible to focus on con- 
straining the heavy quark PDF through measurements 
of this process. 

7.2. Intrinsic Charm 

In the standard global analysis of PDFs, heavy quarks 
do not have free fit parameters associated with them. 
Therefore currently it is assumed that there is no intrin- 
sic charm (IC) or intrinsic bottom (IB) inside the nu- 
cleus. This means that the heavy quark PDF can be 
obtained purely perturbatively through the use of the 
DGLAP evolution equations, having set an appropriate 
initial condition (such as fqix, Qo = mg) - 0). How- 
ever there is no theoretical reason that this should be 
the case and moreover there have been some data that 
point towards the existence of an intrinsic charm com- 
ponent in the proton (Fj large x as measured by EMC 
f2\). As a result there are non-perturbative models pre- 
dicting the size and shape of the intrinsic heavy quark 
component ^ |4l which are used in PDF global fits ||5| 

In Fig. [T|the difference between different IC PDFs 
and the radiatively radiate charm (CTEQ6.6M) is dis- 
played. Both the BHPS (dashed red line) and Meson 
Cloud model (dotted green line) display a similar be- 
havior where they peak at high x, while being very sim- 
ilar to the radiatively generated charm (solid black Une) 



^Currently there are fits available only for IC. The IB/IC ratio is 
expected to be proportional to m^/m^. 
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at lower x. The purely phenomenologically inspired 
sea-like model (dash-dotted blue line) produces a PDF 
which is larger than the radiatively generated one in a 
broad x range. In order to constrain precisely the HQ 
PDF it is important to measure processes sensitive to it 
such as y + 2- How these different models affect the 
y + Q cross-section is presented in the next section. 



the radiatively generated charm and the sea-like charm 
cross-section, however it still undershoots the data at 
very high pj- Therefore in order to confirm or disprove 
the existence of IC more data for processes sensitive to 
IC is needed. In the next subsection the effects of IC on 
the direct photon and heavy quark cross-section, and the 
possibility of constraining IC at the LHC are presented. 



Different Charm PDFs,c(x.Q=40GeV) — CTEQ6.6M 
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Figure 1 : Comparison between the radiatively generated charm PDF 
(CTEQ6.6M, sohd black line), charm PDF including IC: under the 
HHPS model (CTEQ6.6C1, dashed red line), under the Meson Cloud 
model (CTEQ6.5C4. dotted green line), under the sea-like model 
(CTEQ6.6C3, dash-dotted blue line) 



Figure 2: da-y+Q IdpTj-. for charm at NLO (black solid line), at LO 
(black dash-dotted line), for bottom at NLO (red dashed line), at LO 
(red dotted line). 



2. Predictions: hadron-hadron collisions 

2.1. Tevatron 

The differential cross-section vs pjy for y -i- c and 
y + b production at the Tevatron is presented at both 
NLO and LO in Fig. |2] From there we see that the 
difference between dcr'^^'^ Idpjy and dcr^^'^ Idpjy at LO 
stays almost constant, whereas it decreases at NLO with 
growing pry The reason for this is that at the Tevatron, 
due to the abundance of valance quarks and anti-quarks, 
the annihilation subprocess (q + q ^ Q + Q + y) starts 
to dominate the cross-section, and this subprocess is ex- 
actly the same for both y+c and y-i-/? production. When 
one looks at the comparison between theory and exper- 
imental measurements (performed by the D0 collabo- 
ration 16J ) in Fig. [3] the agreement for y + b is very 
good, whereas for y -i- c the data overshoot the theoret- 
ical predictions at high pr- One possibility to try and 
correct for this discrepancy is to utilize IC PDFs in the 
theoretical predictions. In Fig. fflthe data over theory 
ratio is presented, as well as the ratios jg-cTEQeAuf^p^ ; 

^^cTEQiJtfl'Jp^ ■ As expected the sea-like cross-section 
overshoots the data at low pj, but undershoots it at high 
Pt- Conversely the BHPS cross-section, follows the 
data trend, and provides a better description than both 
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Figure 3: Comparison between theory - y + c (solid black curve), y + & 
(dashed red curve) and D0 data -y + c (black circles), y + b (red 
squares). 

2.2. LHC 

At the LHC, already at 7 TeV the x probed is small 
{x ~ ptI VS). Therefore in order to be sensitive to the 
BHPS IC, which shows up at large-jc (Fig. [T]), one has 
to probe forward rapidities in order to be sensitive to 
it. Utilizing experimental cuts appropriate for the CMS 
detector |7, 8| and focusing on the forward rapidity re- 
gion (L566 < \yy\ < 2.5, L5 < lyg] < 2.0) for both 
photons and c-quarks, we present in Fig. |5] the predic- 
tions for do-^^ydpTy at Vs" = 7 TeV. The BHPS cross- 
section (green dashed line) is similar to the radiatively 



2 







- _ H=0.5p.,.^ 




— H=2p-i-^ 




-- BHPS 




■■■■ Sea-like 













_ ly^kl.4442. lyykO.5 


' BHPS/CTEQ6.6M 


_ ly^kl.4442. 1.5<lyyk2 




1.566<lyl<2.5,ly^k0.5 




__ 1.566<ly^l<2.5, 1.5<lyQl<2 







Figure 4: Data over theory ratio for -y + c production (red squares), the 
theoretical predictions are for radiatively generated charm including 
the scale dependence (black lines), BHPS IC (red dashed line), sea- 
Uke IC (blue dotted line). 



Figure 6: Ratio of BHPS IC y + c cross-section over one using ra- 
diatively generated charm for different restrictions on the photon and 
charm rapidity. 



generated one (black solid curve) at low pr, but as pr 
increases (and correspondingly x increases) so does the 
difference between the two curves, thus providing a pos- 
sibility for testing the existence of BHPS IC. To portray 
the need for measurements at forward rapidities the ra- 
tio between the BHPS cross-section and the radiatively 
generated one for several different combinations of the 
photon and heavy quark rapidities is presented in Fig. 
[6] There the ratio at larger pr increases from ~ 1.1 
for \yy\ < 1.4442 & lygl < 0.5, to up to ~ 2.2 for 
1.566 < \yy\ < 2.5 & 1.5 < Ijq] < 2.0. 
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Figure 5: dcf^'^ jdpjy at CMS for forward rapidity for radiatively 
generated charm (black solid line) + scale dependence (black dashed 
curves), BHPS IC (green dash-dotted line), sea-like IC (red dotted 
line). 



3. Predictions: hadron-ion collisions 

It is also interesting to study the IC effects in pA 
collisions. However in order for such a study to be 
performed, the gluon nuclear PDF (nPDF) should be 
as well known as in the free proton. Unfortunately 



this is not the case, and the gluon distribution is only 
very weakly constrained by NMC data on F^ix, Q^) 
and /^2"/^2 2^-* Q large uncertainty in the nu- 
clear gluon as well as the variations in predictions of 
different fits (nCTEQ E [TOl EU, HKN07 (red dash- 
dotted band) US, EPS09 HI) is presented in Fig. |7] 
in the form of the gluon nuclear modification factor 
Rg(x,Q) = g''"'''(x,Q)/g''{x,Q). It is quite clear that 
there is a strong need for data constraining the gluon 
nPDF In the general framework, where there is no IC 
considered, and the charm distribution is solely based 
on the gluon PDF, the charm nuclear modification, Rc, 
follows Rg quite closely as shown in Fig.|8] Since cr''+'^ 
is quite sensitive to the gluon PDF as well as the charm 
PDF as shown in Fig. |9] where the dominating subpro- 
cesses at = 8.8TeV for p - A collisions are gg and 
gQ initiated, it is a very useful process for constrain- 
ing the gluon nPDF. Focusing on the central rapidity re- 
gion in ALICE will isolate any effects caused by the 
possible presence of the theoretically favored BHPS- 
IC, which will appear at large-x and forward rapidity, 
Fig.|6] Therefore we can expect that the nuclear correc- 
tions to 0"''+'^ will follow the ones for Rg and Rc corre- 
spondingly. In Fig. [To] the nuclear modification factor, 
7?"^' ' - T— r-7?r^^^7-^ , to the direct photon and heavy 

Af.i,da-/dpTyipp)' ^ ■> 

quark cross-section is shown. It can be seen clearly 



by comparing Fig. 10 and the boxed region in Fig. [T] 
that the nuclear modification factor to the cross-section, 
R"^ ' , follows closely the nuclear modifications for the 
gluon nPDF (Rg), in the region of x probed at the LHC; 
for more details see ifTSll . Therefore this process is an 
excellent candidate for constraining the gluon nuclear 
distribution. 



*This is true for most nPDF fits, while EPS09 also includes data 
on if' production at RHIC. 
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Figure 7: The gluon nuclear modification factor for nCTEQ (solid 
black line), HKN07 (red dash-dotted band), EPS09 (dashed blue 
band). 



Figure 9: Subprocess contributions to cf*'^ at p-A collisions at AL- 
ICE. 
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Figure 8: Comparison between the charm and gluon nuclear modifi- 
cation factors for nCTEQ (solid black line), HKN07 (red dash-dotted 
band), EPS09 (dashed blue band), the two shaded regions represent 
the X probed at LHC (first band) and at RHIC (second band). 



4. Conclusion 

It was shown that through direct photon production in 
association with a heavy quark we can study and con- 
strain the heavy quark PDFs. More specifically mea- 
surements at the Tevatron of this process have indicated 
a possible existence of IC that seems to favor the BHPS 
model. Further tests at the LHC of this process espe- 
cially at forward rapidities should be able to constrain 
or disfavor the presence of intrinsic charm in the pro- 
ton. While such a constrain inside the nucleus is also 
very important, one has to have a quite precisely deter- 
mined gluon nPDF in order to do so. Since this is not 
the case, we have shown here that direct photon + heavy 
quark production at pA collisions is also a great process 
for constraining the nuclear gluon PDF. 
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Figure 10: Nuclear modification factor for y + c at ALICE, using 
nCTEQ (black solid line), EPS()9 (blue dashed line) + error band, 
HKN07 (red dash-dotted line) + eiTor band. 
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